How to observe a non-Kerr spacetime 
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We present a generic criterion which can be used in gravitational-wave data analysis to distin- 
guish an extreme-mass-ratio inspiral into a Kerr background spacetime from one into a non-Kerr 
background spacetime. The criterion exploits the fact that when an integrable system, such as the 
system that describes geodesic orbits in a Kerr spacetime, is perturbed, the tori in phase space 
which initially corresponded to resonances disintegrate so as to form the so called Birkhoff chains on 
a surface of section, according to the Poincare-Birkhoff theorem. The KAM curves of these islands 
in such a chain share the same ratio of frequencies, even though the frequencies themselves vary 
from one KAM curve to another inside an island. On the other hand, the KAM curves, which do 
not lie in a Birkhoff chain, do not share this characteristic property. Such a temporal constancy of 
the ratio of frequencies during the evolution of the gravitational-wave signal will signal a non-Kerr 
spacetime which could then be further explored. 
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1. INTRODUCTION 

While ground-based gravitational-wave detectors are 
already in operation, and are trying to detect gravita- 
tional waves from stellar mass compact objects, LISA 
-the space-borne detector [l| which is planned to be 
launched during the forthcoming decade- is expected 
to observe much more massive sources of gravitational 
waves with high signal-to-noise ratio. Such signals will 
offer us an opportunity to map the strong field around 
massive astrophysical objects ^\3^, and to test the con- 
ventional wisdom, partly supported by astrophysical ob- 
servations 0, 0|, according to which the massive com- 
pact objects harbored at galactic centers should be highly 
spinning Kerr black holes. Such an astrophysical fact is 
further enforced by the no-hair theorem, which states 
that when a black hole is formed, its multipolc moments, 
and consequently its gravitational field, are determined 
by only two parameters; its mass and spin (assuming its 
charge is negligible). In this paper we suggest a clear and 
generic observable signal that distinguishes an extreme- 
mass-ratio inspiral (EMRI) related to a non Kerr black 
hole spacetime from a corresponding Kerr one. 

Our suggestion comes from the fact that a Kerr space- 
time leads to an integrable system that describes geodesic 
orbits, while any other generic, stationary, and axisym- 
metric non-Kerr spacetime is not expected to be inte- 
grable. Thus assuming that a slightly perturbed Kerr 
metric, that supposedly describes the neighborhood of 
an axisymmetric rotating compact object, is governed by 
a non integrable system of geodesic equations of motion, 
any qualitative new characteristics of the orbit of the test 
body, which could in principle be observed through grav- 
itational waves, can be used to identify such a source. 

According to the KAM theorem Q for dynamical sys- 
tems almost all KAM tori in the phase space of a per- 



turbed integrable system are not destroyed; they simply 
become slightly deformed. However, among the KAM 
tori of the initial integrable system, there are the so 
called resonant tori that are characterized by commensu- 
rate ratios of frequencies. In the perturbed system these 
tori disintegrate, and according to the Poincare-Birkhoff 
theorem j8| they form a chain of islands, on a surface 
of section, inside which the ratio of the corresponding 
frequencies remains equal to the rational number of the 
corresponding initial resonant torus. The width of this 
chain of islands is a monotonically growing function of 
the perturbative parameter that measures the system's 
deviation from the corresponding integrable one, at least 
for small perturbations. 

In our case, an EMRI in a perturbed Kerr spacetime 
will be described by an adiabatically changing geodesic 
orbit, which will sweep a finite range of KAM tori in the 
course of time, while the frequencies of the orbital oscil- 
lations on the polar plane (the plane that passes through 
the axis of symmetry of the massive object and rotates 
along with the low-mass body) will change continuously. 
When the system enters a Birkhoff chain of islands on 
a surface of section, the ratio of frequencies will remain 
constant, while the frequencies will continue to change. 
Therefore the appearance of a plateau in the ratio of fre- 
quencies during the evolution of an EMRI will definitely 
signal the presence of a non-Kerr spacetime. These fre- 
quencies will be encoded in the gravitational wave that is 
radiated from the corresponding source. Thus they can 
in principle be monitored. The question is how probable 
is for an orbit to cross such a chain of islands during its 
evolution that is monitored by the LISA detector. We ar- 
gue [9] that this happens quite often, since the orbits of 
EMRI's that develop in the neighborhood of the super- 
massive compact object at galactic centers are initially 
quite eccentric and inclined [lOj . Depending on the value 
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of its physical parameters, the system will eventually en- 
ter the region of such an island and then the plateau 
in the ratio of frequencies will show up. Furthermore, 
since the rational ratios of frequencies are dense in any 
finite interval, during the whole inspiral phase there is 
a non-negligible probability for observing more than one 
plateaus in frequency ratios. It should be noted though 
that only the low-number resonances are probable to be 
observed since the width of the corresponding islands is 
usually much lower for the higher resonances. 

In order to exhibit our criterion for observing a non- 
Kerr spacetime, we have used a Manko-Novikov (MN) 
metric [U as an example of a non-Kerr metric. The 
specific MN metric is characterized by one parameter q, 
besides the mass and the spin parameters, which mea- 
sures the deviation of its quadrupole moment from the 
corresponding Kerr metric. This particular metric was 
analyzed in [12] and was found to behave like an inte- 
grable system with respect to regularity of its curves in 
the outer allowed region of phase space. Based on the 
general applications of the KAM theorem, we performed 
a more thorough analysis of the orbits in the same met- 
ric and found that the expected Birkhoff chains of islands 
are present on a surface of section, although they are very 
thin to be observed in a coarse study of orbits. 

Also we investigated the effect on the ratio of frequen- 
cies of polar oscillations when an orbit is moving adia- 
batically in and out of such an island while the corre- 
sponding EMRI radiates energy and angular momentum 
away. We have shown that the plateau in the ratios of 
frequencies could last for a few hours to a few weeks, 
depending on the masses involved in the EMRI, as well 
as the value of the q parameter of the background space- 
time. Actually, the aforementioned plateaus will be more 
distinctly observed in lower mass-ratios {fi/M < 10~^ 
where ^, M are the reduced and total mass of the bi- 
nary) since then the system spends more time within an 
island. These effective plateaus could be somehow inte- 
grated in the schemes used in data analysis of LISA to 
look for non-Kerr EMRI's. 

Finally, we should emphasize that the existence of a 
plateau should be a generic result for any perturbed Kerr 
metric, therefore the MN metric which was assumed for 
the background metric is an example that does not imply 
any restrictions in the general application of the proposed 
observable criterion. 



2. RESONANCES IN A NON-KERR MN 
SPACETIME 

The "bumpy" spacetime, which was used by in [T3| to 
explore the orbits in a non-Kerr metric, is a family of 
solutions of the Einstein equations in vacuum, that has 
been built on the foundations of a Kerr metric and thus 
it can be transformed into a pure Kerr metric pT| . For 
the specific MN metric this could actually be done by set- 
ting the q parameter (the quadrupole deviation param- 



eter) equal to zero. The higher mass and mass-current 
multipoles are also affected by q [T^ . 

The Kerr metric is a very special metric for two main 
reasons, (a) It is a highly symmetric metric since apart 
from the three integrals of geodesic motion that all axially 
symmetric, stationary, and asymptotically flat metrics 
share, it is characterized by an extra integral of motion, 
the so called Carter constant, (b) It is a metric related to 
realistic physical objects. According to the "no hair the- 
orem" , all spinning objects that collapse to form a black 
hole are described by such a metric. Since highly com- 
pact objects that move very close to each other are very 
powerful sources of gravitational radiation, an EMRI of a 
stellar-mass compact object around a massive Kerr black 
hole is a highly promising source for the LISA detector. 




FIG. 1: The surface of section of the outer region of allowed 
orbits on the p, p plane for the parameter set E — 0.95, Lz = 
3M, X = 0.9, q = 0.95 of MN (for the definition of these 
quantities see [Ij]). 

In order to look for qualitative new characteristics in 
non-Kerr compact sources that could exist in nature, we 
have computed numerically the geodesic orbits in a MN 
metric. Since this is a metric which is built on the ba- 
sis of a Kerr metric, while it does not share its special 
symmetry that leads to Carter constant, it could be con- 
sidered as a perturbed Kerr; therefore it is not expected 
to be described by an integrable system, in contrast to 
what happens in a Kerr spacetime. The fact that this is 
not an integrable system could already be deduced from 
the fact that the orbits in the inner allowed region of the 
MN spacetime exhibit chaotic behavior, at least in all 
cases studied in [l^]. On the other hand the apparent 
regularity of orbits in the outer region is just an outcome 
of the KAM theorem, as noted in In this region 

(where the gravitational field is weaker) the MN space- 
time looks like a perturbation of Kerr and thus most of 
KAM tori just modify their shape. However, the reso- 
nant tori disintegrate and intersect a surface of section 
forming regions of finite thickness, instead of single closed 
KAM curves. These regions form sets of islands, known 
as Birkhoff chains of islands; each one of them bounded 
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externally and internally by K AM curves [l^ . The exis- 
tence of such chains of islands is characteristic of a system 
that is nearly integrable. A thorough analysis of geodesic 
orbits in the MN spacetime revealed two such Birkhoff 
chains of islands that correspond to the resonances 2/3 
and 1/2. Discovering such islands in a Poincare surface 
of section is not an easy task since these islands are not 
wide enough (cf . Fig. [T|) ; a very fine sweep of initial con- 
ditions is needed to get such an island. Fortunately, there 
is an alternative tool to approach the corresponding or- 
bit much faster, the "rotation number" . This number 
can be easily computed for every orbit and by changing 
the initial conditions and monitoring the corresponding 
rotation number, we could arrive at the desirable rational 
value of the rotation number, by a method analogous to 
a root-finding algorithm through bisection. 

A phase orbit in an integrable system (like the Kerr 
case) is wound around a non-resonant torus filling the 
whole torus densely in the course of time, while for a 
resonant torus the orbit is periodic, repeating itself after 
a few windings f^. Thus, on a surface of section which 
intersects transversally the corresponding tori (we have 
used the plane z = 0, on which we mark the intersecting 
points when i > 0), the crossing points constitute a set 
of points that densely fill a closed curve in the former 
case, and a finite number of fixed points in the latter 
case. This number is the number of oscillations along 
the spatial axis that intersects transversally the surface of 
section (the z-axis in our case) that correspond to a finite 
number of oscillations along the other spatial axis (the p 
axis). When the system deviates from being integrable, 
these fixed points of the periodic orbit evolve to an equal 
number of Birkhoff islands. 

By means of the sequence of crossing points described 
above it is easy to define the aforementioned rotation 
number. Either for an integrable or for a non integrable 
system of 2 degrees of freedom, this is defined as fol- 
lows: Let A be the fixed point on a surface of section 
and Bi be the i-th point of intersection of the phase or- 
bit with the surface of section. The vector ABi rotates 
through an angle A^^ as it moves from the i-th inter- 
secting point to the next one. The rotation number is 

1 1 ^ 

v ~ — lim — A6i (for a review see [131). As ex- 

2tt n^oo N ^ 

i=l 

plained above the rotation number is the fraction of the 
two corresponding fundamental frequencies of the sys- 
tem, each multiplied by an integer. Thus for a geodesic 
orbit in a MN metric it will he i' — ^^^^ with m, n integer 
numbers. Whenever the two frequencies are commensu- 
rate to each other, the rotation number is rational. Actu- 
ally the resonant tori of the integrable system that form 
Birkhoff islands on a surface of section when the system 
gets perturbed, are the ones with commensurate frequen- 
cies, and consequently with rational rotating numbers. 
Thus we can alter the initial conditions until we obtain a 
specific rational value for v and thus reveal the location of 
a Birkhoff island. For a given perturbed system, the ac- 



tual width of the Birkhoff islands is often higher for lower- 
integer ratios of frequencies, while for higher-integer ra- 
tios the islands are so thin, that they are very difficult 
to be discerned. In our case of geodesic orbits in MN we 
have searched for Birkhoff islands when fp/fz = 2, or 1, 
which correspond to = 2/3, or 1/2, respectively. On the 
other hand the thickness of the particular islands grows 
with the perturbation parameter though not linearly, 
in agreement with simple perturbed integrable Hamilto- 
nian toy models that have already been studied in the 
literature [ll,[3[ll. 



3. EVOLVING ORBITS DUE TO RADIATION 
REACTION 

Each chain of Birkhoff islands, separate the surface of 
section in an interior and an exterior region. While one 
moves from the former to the latter region, the rotation 
number, and consequently the ratio of frequencies, drops 
monotonically, except if we land on the thin window of 
an island; then the ratio of frequencies does not change, 
even though the frequencies themselves change. Thus 
by analyzing the spectrum of a gravitational wave signal 
and monitoring the ratio of the fundamental frequencies 
(cf. Fig. [2), whenever a plateau shows up, it will sig- 
nal a non-Kerr type of background. The plateau will be 
more prominent (i) if the deviation from Kerr is larger, 
(ii) if the ratio of mass is more extreme, since then the 
deviation from geodesic orbit, due to energy and angu- 
lar momentum loss, is slower and the system needs more 
time to cross an island, and (iii) the island corresponds to 
a resonance of low-integer ratios, since these resonances 
correspond to wider islands. 

We have run such adiabatic orbits in MN spacetimes 
numerically, following approximate formulae for the evo- 
lution of energy and angular momentum along the axis 
of symmetry (for a generic spacetime there is no reliable 
way to compute the losses due to radiation) [l^ . Our nu- 
merical explorations show that these plateaus are visible 
for ratios of masses below some threshold that depends 
on the resonance we are dealing with and the value of q 
assumed for the MN metric. For example for q = 0.95 
the higher ratio of masses for which we can tell the pres- 
ence of such a plateau is /i/M = 10~*. If the ratio of 
masses is above these thresholds, the system evolves so 
fast that the corresponding plateaus are not discernible. 
Below the threshold, the actual duration of the plateau 
is more extended for lower values of fi/M, although it 
varies a lot depending on the specific trajectory of the 
phase orbit through the Birkhoff island. The aforemen- 
tioned thresholds are actually within the relative range of 
masses expected for an EMRI signal detectable by LISA 

Another crucial point with respect to observability of 
the plateau effect is the following: Since the chains of 
Birkhoff islands are numerous, and the phase orbit is 
forced to move adiabatically from one KAM curve to the 
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FIG. 2: The evolution of the ratio of fundamental frequen- 
cies /p//z for a MN EMRI with q — 0.95 and ratio of masses 
fj./M = 5 X 10~^. The initial conditions were chosen to get 
quite a long plateau. The fundamental frequencies were com- 
puted numerically from the Fourier analysis of the orbital 
motion in time intervals of length 5000M. The oscillations of 
the ratio -mainly at the plateau interval- are artifacts due to 
Fourier analyzing a finite part of the orbit. 

next, it is forced to pass through many such islands. Thus 
while the effect is always present, an actual system will 
exhibit an unambiguous plateau whenever it has suffi- 
cient time to cross a strong resonance, like the 2/3 or the 
1/2 one. This happens if the geometric characteristics 
of the orbit when it enters the window of sensitivity of 



LISA are such that during the evolution of the orbit, it 
will cross one of these resonances. Thus the orbit should 
have initial eccentricity and inclination within a suitable 
range. These ranges are quite wide; therefore it is an- 
ticipated that a large fraction of such suitable EMRI's 
will leave their imprints on their signal through an ap- 
parent plateau of the ratio of the observed frequencies. 
Moreover even if the evolution of the signal is such that 
it correlates well with a Kerr EMRI, we could focus our 
search in the particular period when the ratio of the cor- 
responding peaks in the spectrum that are related to the 
radial and precessional oscillation are close to a resonant 
ratio {fp/fz = 2 for the 2/3 resonance and 1 for the 1/2 
resonance) . A statistically important persistence of such 
a ratio would be a clear "smoking gun" for a non-Kerr 
metric. 

Finally, we should note that a possible positive sig- 
nal from a non-Kerr EMRI, could be further explored 
for other consequences of such a peculiar source, like the 
instabilities that are expected to show up before the cor- 
responding final plunge [l2|, or a possible transit of the 
orbit to chaotic behavior through entrance in the interior 
region of allowed orbits of a MN-likc metric |9|, 112] ■ 
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